Fluctuation induced non-canonical BCS states: A mechanism for 
pseudogap 
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Abstract. - We pose the question of what effect the statistical fluctuation causes if it induces 
a non-unitary evolution. We apply this idea to the BCS model and study fluctuation around 
the mean-field average. We find that, dynamics of the thermalization influences the equilibrium 
besides the non-unitary evolution, and the resulting equilibrium state is no longer the canonical 
one. The pseudogap phenomenon can exist in this model. 
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Introduction. — Consider equilibrium states of a sys- 
tem governed by master equation 



dp 
dt 



-i[H,p]-7iA,[A,p]l 



(1) 



where A is an operator of the system. The equation is a 
simplified Lindblad equation |T] with the operators speci- 
fied as a Hermitian one. If an energy eigenstate \E){E\ 
is not an eigenstate of ^1 , it is in general not time- 
independent, d{\E){E\)/dt 7^ 0. In this case, the canon- 
ical state p cx qt^/'^bT equilibrium state, 
and a non-canonical one is required. Suppose the goal 
of thermal relaxation is the canonical equilibrium ensem- 
The thermalization ceases when the goal is achieved. 
When the system is in a non-canonical equilibrium state, 
however, the thermalization still works. So the dynamics 
, of thermalization is also a factor in determination of the 
equilibrium. 

The BCS theory is a non-interaction quasiparticle 
model. The subadditivity [5] of entropy for a composite 
system p^^ indicates S{p^^) < S'(TrA[p'*^] «) Tt:b[p^^])- 
When there is no interaction between the two subsystems, 
this inequality can be generalized to the free energy as 
F{p^^) ^ F{Tya[p^^] ® Trsb-^^]). As long as A does 
not mix states of different momentum, according to the 
principle of minimum free energy, a BCS state after ther- 
mal relaxation must take the form p — ®kpk [3^, where 
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^The argument is given in the Appendix. 



Pk is a density matrix in the subspace of momentum k. 
So we can handle each subspace separately, which allow 
us to study the energy variation caused by a gap for each 
momentum. We assume that a gap emerges only in the 
subspaces where the energy is reduced by the gapping. For 
the gapped momentums, we adopt the Cooper approxima- 
tion, a constant Afc. 

Here, a subspace of momentum k means the one 
spanned by the ground state component, two excitons and 
the double exciton: 
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where | ) is the bare vacuum, b'^^,b'^f,^ the Bogoliubov 
creation operators and c^k^,cl-^ creation operators of an 
electron. The subspace spanned by states of fc | and —fet 
are referred to by —k. 

To have a sense of how we can deal with a subspace 
separately, let us derive the Fermi-Dirac distribution for 
the BCS model according to the principle of minimum 
free energy. A density operator in a subspace can be 

expressed as pk = Smn "™nl"^) ("-I "^ith "mn = "nm- 
Since the time-independent states of von Neumman equa- 
tion are energy eigenstates, an equilibrium state has form 



Pk 



E 



n=0 ■ 



|n)(n|. We assign variables as aoo = 
p_,a33 = p+,aii = a22 = (1 - - p+)/2 and set 
the energy of the ground states zero. Then we have 



p-1 
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= = E^/2 = Ek = Vt^2 ^ A2,w = Efc - ep- By 
minimizing the free energy, we obtain the probabilities 

p^ = {e^f'^- + l)-\ (3) 

where (3 = ksT. The Fermi-Dirac distribution is given by 

(l-p_+p+)/2 = (e'3^^'= + l)-i 

Non-unitary evolution and thermalization. A 

many-body Hamiltonian can not lead to non-unitary evo- 
lution by itself, because no environment is involved. For a 
part of this many-particle system, however, the other parts 
actually play a roll as an environment with fluctuations, 
which can cause decoherence d]. Consider the subspace 
of momentum fe, Hamiltonian for which is given by 

Hk = Ek{b+^bkt + btkib-ki) + P, (4) 
where P is the fluctuation around the mean field average 



P 



4t'=ifci("3) Xl'^'^-'^'^^'^'t " (c-fc'iCfe't>) + c-c. (5) 



From C-kiCk-r = (Afe/£'fe)(l - n^ki - "-fct)' P actually 
a fluctuation concerning the particle number in gapped 
states, and it sets off when a gap emerges. To have an 
explicit expression of the non-unitary dynamics cause by 
this fluctuation, we assume the fluctuation can be ap- 
proximated by a Gaussian noise. The trivial equality 
(c-fc'4,Cfc'-f-) — (c_fe'iCfc't) — indicates that the fluctua- 
tion is a non-biased noise. We specify the phase factor so 
that Afc is a real number. Then the subspace is coupled 
to the noise though operator A = c^^c'^/^^ + c^kiCk^ with 
coupling —g. In the "interaction picture" the non- unitary 
evolution has a Lindblad form ^ 



dpi 

dt 



-g'T,[A, [A, pi]] 



(6) 



where Tc = f{t)f{0)dt is the autocorrelation for t — 0, 
at which the gap is opened up. In the Schrodinger picture, 
Eq. dH) takes the form of Eq. 

Since eigenstates of A are different from the energy 
eigenstates, we need to construct a thermalization model. 
Effect of the thermalization is reduction of the free en- 
ergy. To achieve the canonical ensemble, we may expect 
a phenomenological description of the thermalization as 
the following: for off-diagonal elements amn,™ 7^ n, the 
time-dependence obeys 



da„ 



dt 



fmn {Pk^^mn ; 



while diagonal elements a„„ vary with time as 



dttr, 



dt 



-hnn{Pk) 



dF 

ddn.n. 



(7) 



(8) 



Here, fmnipk) , hnn{pk) are positive functions of pfc repre- 
senting thermalization properties of the system, and 



F 



1 

E" + TkBttnn lna„ 



(9) 



is the free energy without considering the off diagonal ele- 
ments. For simplicity, we consider the functions two con- 
stants, 2/ and 2h respectively. The factor 2 is included to 
simplify latter notation. 

The overall dynamics of our BCS model consists of three 
components, that is, the von Neumman equation, the non- 
unitary evolution and the thermal relaxation (_R), i.e. 

^ = -i[H,Pk]-g^T,[A,[A,pk]]+R (10) 

Without the non-unitary evolution, the off-diagonal ele- 
ments vary as 



.it) 



lito] 



(11) 



which implies that they approach zero through relaxation. 
Then the von Neumman term gives zero because energy 
eigenstates are commutable with the Hamiltonian. The 
master equation (fTOl) reduces to Eq. (jH]), where the time- 
independence is just the condition for minimum free en- 
ergy. The solution is the canonical equilibrium state. 



The pseudogap. 

C-ki^k-f acts as: 

M 

Ac 



The operator A 
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= Ac 



> =0. 



(12) 



Action of [A, [A, p]] annihilate the state elements except 
]m){n]{m, n = 0, 3). The non-unitary evolution impose no 
constraints on matrix elements am„(m, n 7^ 0,3). Those 
elements approach zero through thermal relaxation except 
flu and 022. We only need to consider the following matrix 





Pk = 



( floo 



V 030 





an 







0-22 





ao3 \ 



033 / 



(13) 



and solve the equations given by coefficients of the six 
state elements. Among the six equations, the two from 
|1)(1|, |2)(2| have the same solution, and the two from 
|0)(3|, |3)(0| are equivalent because of the Hermitian sym- 
metry. So only four equations are independent. 

It is more convenient to obtain the equilibrium 
state through equations from the coefficients of 
l>( 



which are given by 

\UkVk{am - ass) - {h 
Viul ~ vl){aoo - 033) + 
-iufcWfc(aoo - 033) + (h 
V{ul- vl){aoa- 033) + 



I )( \c-kiCkt, I )( \c-kiCkt and c+^^| )( ]c-ki, 



2r])ukVk{ao3 + 030) 



dF 



OF 



^ daoo " da33 
'2.T])ukVk{ao3 + 030 
, dF 



= 0, (14) 



; dF 
'daoo 



'■da 



33 



0, (15) 



-iEulao3 - iEvla3o + 77(003 - 030) 



UkVk{ 



dF dF 



da, 



00 



da33 
dF 



) = 0, 



dan 



= 0, 



(16) 
(17) 
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Non-canonical BCS states 



where rj — g^Tc- 

From summation of Eq. 

dF 



Tij) and Eq. 



T5l). we have 



■00 



da. 



33 



(18) 



Thus Eq. ([T7| holds according to an = (1 — aoo — a33)/2 
and the chain rule of derivative. Owing to {uk-x) + 
{n-k\) = 1/2— (aoo — 033)72, what concerns our discussion 
is the value of aoo — 0,33- Eq- dl) gives relation 



(19) 



1 + (aoo - a33)^ - 2(aoo + 033) = 0. 
Substituting it in dF/daoo, we have 



dF TP , rr, , 1 + aoo - 033 

-5 = -Ek + TkB In r (20) 

oaoo 1 - (aoo - 033 ) 

We see that Eq. ([14]) is a transcendental equation of 
ooo — «33- To obtain an approximate solution, we con- 
sider the first order expansion of ooo — 033 with limit of 
strong interaction, a large 77, which leads to an approxi- 
mate solution 



floo — 033 — Ek{ 



hfEi + 2r//iA2 



2kBT)- 



(21) 



Our approximation requires aoo — 0,33 ^ 1- From the solu- 
tion, this condition is ensured when A is a small number, 
A ~ 0. When a; ^ 0, a small ^ makes aoo — ^33 
a small number. For cj 9^ 0, a large 77 results in a small 
aoo — 033 ■ It is necessary to note that the solution not only 
applies to the gapped momentum, but also to the states 
ungapped. This is because the non-unitary evolution af- 
fects on all momentums as long as some momentums are 
gapped. For the gapped momentums, A equals the gap 
constant, while A = for those ungapped momentums. 

According to our assumption, the gapping occurs only 
in the subspace where a gap reduces the energy. For an 
infinitesimal gap, the gapped momentum should satisfy 



A2=0 



< 0. With our approximation, the derivative 



is given by 



d{5Ek) 



d(A2) 



A2=0 



1 

2w 



2h 
J' 



(22) 



For momentums near the Fermi surface, we have ~ 0, 
so the gap for them can not be opened. With proper 
parameters / and h, momentums having a distance to 
the Fermi surface can be gapped. Because the ungapped 
states can be excited without overcoming an energy gap, 
the overall state is not in the superconducting phase. Thus 
the pseudogap phenomenon emerges. 

Conclusion. — We studied the non-unitary evolution 
caused by the fluctuation around the mean field average. 
We showed that such an evolution can break the supercon- 
ducting gap for momentums near the Fermi surface and 
lead to the pseudogap phenomenon. Because the equilib- 
rium state is a non-canonical one, dynamics of the thermal 



relaxation also influence its form. In the non-canonical 
equilibrium states, the particle number does not obey the 
Fermi-Dirac distribution even in non-interaction Fcrmion 
systems. 



The author is grateful to X. Sun, C.Q. Wu and Y.L. Ma 
for useful discussions. The work was partially supported 
by the Doctoral Foundation Program of the Chinese Min- 
istry of Education. 



REFERENCES 

[1] LiNDBLAD G., Commun. Math. Phys., 48 (1976) 119. 

[2] Peres A., Quantum Theory: Concepts and Methods 

(Kluwer Academic, New York) 2002, p. 270-275. 
[3] Gu L., Phys. Scr., 48 (2010) 055004. 
[4] ZUREK W. H., Rev. Mod. Phys., 75 (2003) 715. 
[5] Matsuzaki Y., Saito S., Kakuyanagi K. and Semba 

K., Phys. Rev. B, 82 (2010) 180518(R). 

Appendix: Goal of thermalization. Suppose an 
equilibrium ensemble of an N dimensional system is con- 
stituted by |n) = Cnrn\Em), whcrc \Em) (m = 1, 2, • • • , N) 
is energy eigenstates. According to the principle of mini- 
mum free energy, the equilibrium state is given by 



P 



exp(-/3 J2r 



'Era)\n){n\ 



Tr[exp(-/3E,: 



'Em)\n){n\ 



(23) 



Then, taking the partial derivatives with respect to Cnm, 
one can find the canonical equilibrium state (a„„i — 6nm) 
is a local minimum of the free energy, while other states 
are not in general. Of course, there may be other local 
minimums for delicate assignment of the parameters, and 
a rigid argument should show the canonical ensemble is 
the global minimum. In view of the common validity of 
the canonical ensemble, however, we assume it is the goal 
of thermal relaxation. 
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